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We search for the semi-inclusive process using 2.8 fb ^ oi pp coUisions at = 

1.96 TeV recorded by the DO detector operating at the Ferinilab Tevatron Collider. We observe 
26.6 ± 8.4 signal events with a significance above background of 3.2 standard deviations yielding a 
branching ratio of B{B° -> Di*^Di*'') = 0.035 ± O.OlO(stat) ± O.Oll(syst). Under certain theoretical 
assumptions, these double-charm final states saturate CP-even eigenstates in the B° decays resulting 
in a width difference of ArJ^/r^ = 0.072 ± 0.021(stat) ± 0.022(syst). 

PACS numbers; 13.25.Hw, 12.15.Ff, 11.30.Er, 14.40.Nd 



The phenomenon of CP violation is believed to be in- timately tied to explaining the matter dominance in the 
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present day universe [if. CP violation is expected to 
occur in the evolution of neutral particles that can mix 
between different eigenbases. For the system, the 
flavor eigenstates can be decomposed into heavy (H) 
and light (L) states based on mass or into even and 
odd states based on CP. The width differences between 
these eigenstates are defined by AFs = — Th and 
ArJP = r^™" - T°/'^, respectively. These two quantities 
are connected with the possible presence of new physics 
(NP) by AF^ = AF^P cosc/)^, where c/^s is the CP violat- 
ing mixing phase which constrains models of NP. 

In the standard model (SM) a mixing parameter, F12, 
determining the size of the width difference between CP 
eigenstates stems from the decays into final states com- 
mon to both B and B. Since this quantity is dominated 
by CKM-favored tree-level decays, it is practically insen- 
sitive to NP. Due to the hierarchy of the quark mixing 
matrix 0], the width difference is governed by the par- 
tial widths of B'^ decays into final CP eigenstates through 
the b CCS quark- level transition, such as D^D^ 
or i?g J/ipcj). Topologically, the former type of de- 
cay mode is a color-allowed spectator, while the latter 
type is suppressed by the effective color factor. Thus, 
the semi-inclusive decay modes — > ui*^ Di*\ where 
ui*^ denotes either Df or Df*, are interesting because 
they give the largest contribution to the difference be- 
tween the widths of the heavy and light states. The 
other decay modes are estimated to contribute less than 
0.01 to the projected ~ 0.15 value of AF^/Fs d, where 

r,(=l/T,) = (FL+Fjy)/2. 

In the Shifman-Voloshin (SV) limit [4], given by mf, — 
2mc — !■ with Nc — > 00 (where Nc is the number of 
colors), APCP is saturated by T{B° Di*^Di*^). Then 
the width difference can be related to the branching ratio 
of B2 mesons to this inclusive double-charm final state 
by BM 
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where x/ is the fraction of the CP-odd component of the 
decay, F°'i'^/F°™" = x//(l - Xf). Therefore, given the 
CP structure of the final state, AF^^ can be measured 
using the information from branching ratios without life- 
time fits. The irreducible theoretical uncertainty of this 
approach stems from the omission of CKM-suppressed 
decays through the b — > uus transition which is of order 
2\VubVus/V,bV,s\^i-5%. 

In this Letter we report the first evidence for the decay 
B° Di*^Di*\ The study uses a data sample of pp col- 
lisions at y/s = 1.96 TeV corresponding to an integrated 
luminosity of 2.8 fb^^ recorded by the DO detector op- 
erating at the Fermilab Tevatron Collider during 2002 - 
2007. This supersedes our previous study of the same 



final state based on 1.3 fb ^ Q. A similar study based 
on events containing two (j) mesons has been reported b" 
the ALEPH collaboration at the CERN LEP Collider ^ 

This analysis considers the B'^ decay into two Di*^ 
mesons. No attempt is made to identify the photon or 
tt" emanating from the D* decay. We search for one 
hadronic decay to (t>TT and one semileptonic decay 
to where both cj) mesons decay to K~^K~. The 

branching fraction is extracted by normalizing the B'^ 



Di*''Di*'' decay to the B° 



oi*^ decay. 



DO is a general purpose detector [!| consisting of a cen- 
tral tracking system, uranium/liquid-argon calorimeters, 
and an iron toroid muon spectrometer. The central track- 
ing system allows charged particles to be reconstructed. 
This system is composed of a silicon microstrip tracker 
(SMT) and a central fiber tracker (CFT) embedded in 
a 2 T solenoidal magnetic field. Muons are identified 
and reconstructed with a magnetic spectrometer located 
outside of the calorimeter. The spectrometer contains 
magnetized iron toroids and three super-layers of propor- 
tional drift tubes along with scintillation trigger counters. 
Information from the muon and tracking systems is used 
to form muon triggers. For the events used by this anal- 
ysis, the muon from the semileptonic Ds decay satisfies 
the inclusive single- muon triggers. 

Muons are identified by requiring segments recon- 
structed in at least two out of the three super-layers in 
the muon system and associated with a trajectory re- 
constructed with hits in both the SMT and the CFT. 
We select muon candidates with transverse momentum 
Pt > 2.0 GeV/c and total momentum p > 3.0 GeV/c. 

(j) mesons are formed from two opposite sign charged 
particles with px > 0.7 GeV/c in the event assuming a 
kaon mass hypothesis. We require at least one kaon to 
have an impact parameter clearly separated from the pp 
interaction point (primary vertex) with at a minimum 4 
standard deviations significance. The two-kaon systems 
satisfying pt{KK) > 2.0 GeV/c and 1.010 < m{KK) < 
1.030 GeV/c^ are selected as (p candidates. 

The hadronic meson is reconstructed by combining 
the 4> candidate with a third track with px > 0.5 GeV/c 
which is assigned the pion mass. The pion is required 
to have charge opposite to that of the muon. The three 
particles must form a well reconstructed vertex displaced 
from the primary vertex jioj . We require the cosine of the 
angle between the Ds momentum and the direction from 
the primary vertex to the Dg vertex to be greater than 
0.9. For the signal decay chain of a pseudoscalar to a vec- 
tor plus pseudoscalar, followed by the decay of the vector 
to two pseudoscalars, cos is distributed quadratically, 
where 9^ is the decay angle of a kaon in the 4> rest frame 
with respect to the direction of the Ds meson, and hence 
a constraint | cos 0^1 > 0.3 is imposed. 

The i3g Di*^ decay vertex is reconstructed based 
on the momentum and direction of the reconstructed 
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FIG. 1: Invariant mass distribution of the 
(j>Tv system for the Oi*^ fj,u sample. 

The two peaks correspond to the can- 
didates (lower masses) and Ds candidates 
(higher masses). 



FIG. 2; Projections of the two-dimensional maximum likelihood fit onto invariant 
mass spectra of the (a) (jf)7r system from hadronic Ds decays and (b) KK system 
from semileptonic Ds decays. The peaks in both distributions are explored to 
search for the correlation between the two systems. 



hadronic Ds candidate and its intersection with the track 
of an oppositely charged muon. This vertex is required 
to be located between the primary vertex and the Dg ver- 
tex, whereby the individual Bs and Dg vertex displace- 
ments are consistent with a. pp ^ Bg ^ Dg decay chain. 
The invariant mass of the B^ candidate is required to 
be less than 5.2 GeV/c^. We require the daughter parti- 
cles of the i?j meson to be well isolated from other tracks. 
Background is further suppressed using a likelihood ratio 
technique [llj that combines information from the invari- 
ant masses and momenta of the reconstructed particles, 
vertex quality, and the (p helicity angle. 

The (/)7r invariant mass distribution for B^ — > D^*^ 
candidates is shown in Fig. [T] Maxima corresponding 
to the Dg (f>n decay and the <j)TT decay are 

clearly observed. The Dg signal originates from ~ 90% 
semileptonic B^ decays and ^ 10% decays of the type 
B ^ DgD followed by semileptonic D decay. These frac- 
tions are determined from Monte Carlo (MC) simulation 
using the known or estimated branching fractions from 
the PDG [[12] or evtgen [l^. Approximately 2% of the 
events are due to direct charm production pp DD, 
determined by using full simulation and reconstruction 
of DD* candidates. The overall sample composition is 
verified u sing studies of the B lifetime and mixing pa- 
rameters fl4[l5|. 

For the second (p candidate, we search for an additional 
pair of oppositely charged particles in the event imposing 
the same criteria as for the first meson. The two kaon 
tracks are combined with the muon track to produce a 
common vertex for the semileptonic Dg candidate. We 
require the Dg candidate to originate from a common 
vertex to the hadronic Dg candidate to complete the 
Bg Di*'' D^*^ decay. This approach is justified since 
the average transverse decay length of the Dg meson rel- 
ative to the Bg meson decay vertex is ~ 1.0 mm with an 
uncertainty of ~ 0.6 mm. By applying the same selec- 
tion criteria as in the normalization Bg D^*^ /lu decay 
sample, many detector related systematic effects cancel. 



The total invariant mass is required to lie between 4.30 
and 5.20 GeV/c^. 

Correlated production of this double-charm decay, 
where both Dg mesons originate from the same parent 
Bg meson, is then determined by examining the two- 
dimensional distribution of m((/)7r) from hadronic Dg can- 
didates versus m{KK) from semileptonic candidates. We 
perform a maximum likelihood fit to this distribution 
with four components: the correlated DgDg component 
is modeled as the product of signal terms in both di- 
mensions, the uncorrelated components are modeled as 
the product of the signal term in one dimension and the 
background term in the other dimension, and the back- 
ground correlation is modeled as the combination of the 
background terms in both dimensions. Signal and back- 
ground models are expected to be identical with those 
for the Bg —> Di*^ fiu sample, from which the param- 
eters of the signal models are determined. Projections 
of the two-dimensional likelihood fit onto both axes are 
displayed in Fig. ^ The fit returns a yield of 31.0 ± 9.4 
correlated events. 

Three possible sources of background are considered 
in the correlated sample. Direct charm production from 
pp is estimated based on the fraction of prompt charm 
measured directly in the inclusive Dg /ii/ sample, (10. 3± 
2.5)%, along with the decay fraction of the second charm 
quark to a Dg meson and the reconstruction efficiency for 
this decay. Due to a shorter decay length of the charm 
decay, the lifetime requirement reduces its contribution 
significantly leading to an estimate of (1.9 ± 0.5)%. 

The second background source arises from the semilep- 
tonic Bg d'^s'^ (p^LV decay. This can be extracted 
by studying the m{(j3fj,) spectrum. In this variable, 
Bg Dg Dg events tend towards lower values, while 
Bg Di*^(j>^iy events tend towards higher values. 

The third source consists of D^g^D^g^KX 
events. This background can be extracted by study- 
ing the visible mass of all reconstructed daughter par- 
ticles, m{Dg(j)p). The mass tends to have higher values 
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Di*^Di*^ than for B^'° 



for B° 

These backgrounds are estimated with MC samples 
by repeating the fit in three separate regions chosen so 
that mainly one source contributes to each region in the 
m{4>fj.) — m{Ds(j)fJ.) plane. The separate components, the 
signal and the two latter backgrounds, are then extracted 
based on the expected distribution over the three re- 
gions of the three components. We find a signal yield of 
26. 6±8. 4 events originating from the 5° Di*-'!)^*-' pro- 
cess after subtracting the correlated background events. 

The signal is normalized to the total — > ui*^ iiu 
yield taking into account the composition of the sample 
as discussed earlier. The reconstruction efficiency ratio 
between the two samples is estimated from MC to be 
0.082 ± 0.015. This small value results from the softer 
muon momentum spectrum in charm decays as compared 
to bottom decays. The systematic uncertainty in the ra- 
tio contains uncertainties from the modeling of the -B^ 
momentum spectrum, the decay form factors and sam- 
ple composition, and the trigger and reconstruction ef- 
ficiencies. Our efficiency model is verified by compar- 
ing the expected and measured Ds yield and the relative 
50 _^ Di*'^Di*'> to B° Di*'>fiiy yields as a function of 
muon pa- 
using all the above inputs, the branching ratio is mea- 
sured as 

B{B° 

= 0.035 ± O.OlO(stat) ± 0.008(exp syst) ± 0.007(ext), 

where the "ext" uncertainty arises from the external in- 
put branching ratios taken from the PDG ■ This un- 
certainty contributes ^ 45% to the total systematic un- 
certainty (exp syst0ext), which leaves room for further 
improvements in the result. The experimental systematic 
uncertainty accounts for the rest of the total systematic 
uncertainty, containing a 37% component from the re- 
construction efhciency ratio, 11% from the background 
estimation, and 4% from the fitting procedure. All other 
uncertainties are < 1%. 

The probability that the total background would fluc- 
tuate to the measured event yield or higher is evaluated 
to be 1.2 X 10~^ through pseudo-experiments including 
systematic uncertainties. This corresponds to a signifi- 
cance of 3.2 standard deviations. 

Information on the mixing-induced CP asymmetry in 
the B'^ system can be extracted from the branching frac- 
tion measurement through Eq. Since the CP struc- 
ture of the decay is presently not accessible either in 
theory or experiment, several scenarios for different Xf 
values can be considered. In the heavy quark hypothe- 
sis 3 along with the SV limit, the CP-odd component 
of the decay vanishes, leaving the inclusive final state to 
be CP-even, i.e. Xf — 0, with a theoretical uncertainty 
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FIG. 3: Constraints in the ATs — 4>s plane. The solid line rep- 
resents our measurement under the theoretical assumptions 
stated in the text and with xf — 0. Two pairs of lines are 68% 
(dashed) and 90% (dotted) C.L. intervals of AFs for a given 
assumed value of (ps- Contours from the — » J/ip<f> decay 
are the equivalent C.L. regions of (AFs, 4>s) when measuring 
simultaneously both parameters. No theoretical uncertainties 
are reflected in the plot. The SM prediction is represented by 
the thick vertical line. 



measurement assuming the relation AF^ = AF^^ cos ips ■ 
Confidence-level (C.L.) contours from the flavor-tagged 
decay i?g J/iptj) at DO [17] are superimposed. We take 
the mean lifetime of meson from Ref. [12]. 

Furthermore, within the SM framework, the mass 
eigenstates coincide with the CP eigenstates and the ex- 
pression used in the previous studies 0, S] is recovered. 
Our measurement gives 

(*) n(*h 



APCP 2^(5° ^ Di*>Dr) 



of ~ 5% [l6|. This scenario is illustrated in Fig. [31 pre- 
senting the constraint in the AFg — (j>s plane from this 



= 0.072 ± 0.021(stat) ± 0.022(syst). 

This result is consistent with the SM prediction [3] as 
well as with the current world average value [l6| . There- 
fore, if the CP structure of the final state can be disen- 
tangled and the theoretical errors can be controlled, this 
approach can provide a powerful constraint on mixing 
and CP violation in the B'^ system. 

In summary, we performed a study of B^ decays into 
the semi-inclusive double-charm final state using an in- 
tegrated luminosity of 2.8 fb~^ at the DO experiment. 
We see evidence of this process and measure the branch- 
ing ratio as ^(5° ^ Di*^Di*'^) = 0.035 ± O.OlO(stat) ± 
0.011 (syst). Based on this measurement and under cer- 
tain theoretical assumptions, mixing and CP violation 
information in the B'^ meson system are extracted. This 
is the first single measurement that demonstrates a non 
zero width difference in the B'^ system at greater than 
3(7 significance. In particular, in the absence of NP, 
the fractional width difference is derived as AF^^/Fs = 
0.072 ± 0.021(stat) ± 0.022(syst). 



7 



We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAFESF and FUNDUNESF (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); STFC (United 
Kingdom); MSMT and GACR (Czech Republic); CRC 
Program, CFI, NSERC and WestGrid Project (Canada); 
BMBF and DFG (Germany); SFI (Ireland); The Swedish 
Research Council (Sweden); CAS and CNSF (China); 
and the Alexander von Humboldt Foundation (Ger- 
many) . 



[b] 
[c] 
[d] 



[f] 



[h] 
fil 



Visitor from Augustana College, Sioux Falls, SD, USA. 
Visitor from Rutgers University, Piscataway, NJ, USA. 
Visitor from The University of Liverpool, Liverpool, UK. 
Visitor from II. Physikalisches Institut, Georg- August- 
University, Gottingen, Germany. 

Visitor from Centro de Investigacion en Computacion - 
IPN, Mexico City, Mexico. 

Visitor from ECFM, Universidad Autonoma de Sinaloa, 
Culiacan, Mexico. 

Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

Visitor from Universitat Bern, Bern, Switzerland. 
Visitor from Universitat Zurich, Ziirich, Switzerland. 



[fl Deceased. 

[1] A. Sakharov, JTEP Lett. 5, 24 (1967). 
[2] G.C. Branco and L. Lavoura, Phys. Rev. D 38, 2295 
(1988). 

[3] R. Aleksan et al, Phys. Lett. B 316, 567 (1993). 

[4] M.A. Shifman and M.B. Voloshin, Sov. J. Nucl. Phys. 

47, 511 (1988). 

[5] K. Anikeev et al. , |arXiv:hep-ph/020107T] (2002) , Chap. 8. 
[6] I. Dunietz, R. Fleischer and U. Nierste, Phys. Rev. D 63, 

114015 (2001). 

[7] V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

99, 241801 (2007). 
[8] R. Barate et al. (ALEPH Collaboration), Phys. Lett. B 

486, 286 (2000). 
[9] V.M. Abazov et al. (DO Collaboration), Nucl. lustrum. 

Methods Phys. Res. A 565, 463 (2006). 
[10] J. Abdallah et al. (DELPHI Collaboration), Eur. Phys. 

J. C 32, 185 (2004). 
[11] G. Borisov, Nucl. Instrum. Methods Phys. Res. A 417, 

384 (1998). 

[12] C. Amsler et al, Phys. Lett. B 667, 1 (2008). 

[13] D.J. Lange, Nucl. Instrum. Methods A 462, 152 (2001). 

[14] V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

97, 241801 (2006). 
[15] V.M. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

97, 021802 (2006). 
[16] E. Barberio et al. (Heavy Flavor Averaging Group), 

lajXiv:0808.1297l [hep-ex] (2008). 
[17] V.M. Abazov et al. (DO Collaboration), Phys. Rev. 

Lett. 101, 241801 (2008). This measurement establishes 

APs > with a significance of 2.4 standard deviations. 
[18] A. Lenz and U. Nierste, JHEP 0706, 072 (2007). 



